AD—A096  194 

UNCLASSIFIED 

I . 

I 


ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
DIFFUSION  EQUATION  SOLUTION  SEQUENCES. (U> 

JAN  81  J  F  POLK 
ARBRL-TR-02281 


ABERD— ETC 


NL 


SBXE-AD-E430  572 


/L 


TECHNICAL  REPORT  ARBRL-TR-02281 

IFFUS ION  EQUATION  SOLUTION  SEQUENCES 

John  F.  Polk 

January  1981 

DTIC 

^ELECTE| 

% 

MAR  1  1  1981 

w 

B 

US  ARMY  ARMAMENT  RESEARCH  ANO  DEVELOPMENT  COMMAND 
BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


toorovad  for  public  rtlotsoi  distribution  uni  Ini  tad. 


Destroy  this  report  when  it  is  no  longsr  needed. 
Do  not  return  it  to  the  originetor. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  aay  be  obtained 
from  the  National  Technical  Information  Service, 

U.S.  Department  of  Cowmtcb  ,  Springfield,  Virginia 
22151. 


The  findings  in  this  report  are  not  to  be  construed  as 
sa  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 

Tht  um  of  trad $  name*  or  mtmufaaturort '  merit  in  tki*  import 
dot*  rot  oonstitut*  irdoroomtnt  of  mg  oommroial  produet. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (When  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 

TECHNICAL  REPORT  ARBRL-TR- 02281 


4.  TITLE  (and  Submit) 

DIFFUSION  EQUATION  SOLUTION  SEQUENCES 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR  (•) 

John  F.  Polk 


8.  CONTRACT  OR  GRANT  NUMBER/*) 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

USA  Ballistic  Research  Laboratory 
ATTN:  DRDAR-BLT 

Aberdeen  Proving  Ground,  Maryland  21005 


II.  CONTROLLING  OFFICE  NAME  AND  ADORESS 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


RDT!,H  ll.Ktl  102AH43 


tz.  REPORT  DATE 

JANUARY  1981 


■ 

I 


USA  Armament  Research  and  Development  Command 

USA  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BL  Aberdeen  Pr.  Gds.,  MD  21005 


.  MONITORING  AGENCY  NAME  A  ADDRESS///  dUftrtnt  Iron  Controlling  Olllet)  IS.  SECURITY  CLASS,  (oi  tbit  report) 

Unclassified 


13.  NUMBER  OF  PAGES 

49 


15a.  DECLASSI  FI  CATION /DOWN  GRADING 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  (ol  thlt  /Import) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Ihm  ebetrect  an  farad  In  Block  JO,  II  dlllerent  horn  Report) 


If,  KCY  WORDS  (Contlnum  on  rororao  mldm  II  ntcvtiary  and  Idontlty  by  block  numbor) 


Parabolic  Differential  F.quations 
Diffusion  Equation 
Cauchy  Problem 
Boundary  Value  Problem 
Asymptotic  Expansions 


ao.  ABSTRACT  (Centmue  as  reeerm 


Short  Time  Behavior 
Discontinous  Data 
Boundary  Layer 
Interior  Layer 
Stretched  Variable 


r  Identity  by  block  nuaikarj  (  j  f  p  1 


A  methodology  is  described  for  deriving  expansions  for  the  solutions  of  parabolic 
differential  equations  which  are  asymptotically  valid  for  small  times.  Both  the 
Cauchy  problem  and  the  first  initial-boundary  value  problem  are  studied  and  par¬ 
ticular  attention  is  given  to  cases  in  which  the  data  is  discontinuous  or  incom¬ 
patible  with  the  equation.  The  basic  procedure  is  called  the  DESS  (Diffusion 
Equation  Solution  Sequence)  method  and  all  expansions  are  obtained  in  explicit 
form  using  certain  special  functions  which  arc  briefly  discussed.  Error  estimates 
are  rigorously  stated  but  their  proofs  are  not  included  in  the  nresc 


EDITION  OF  1  NOV  •»  IS  OBSOLETE 


UNCLASSIFIED 

SE CUR1TY  CLASSIFICATION  OF  THIS  PAGE  (When  Dote  Entered) 


■  Item.  ***.VK**4MFT4 „ 


TABLE  OF  CONTENTS 


«*  •  *.  « 


tmrnrnimm/mmm 


l.  INTRODUCTION  AND  PRELIMINARIES 

In  this  report  we  study  the  parabolic  differential  equation 
D.u  =  a(x)D2  u  +  b(x)D  u  +  c(x)u  (1.1) 

t  X  X 

and  describe  a  general  method  for  constructing  asymptotic  ex¬ 
pansions  for  its  solutions,  accurate  for  small  times;  both  the 

Cauchy  problem  and  the  first  initial-boundary  value  problem 
will  be  considered.  The  approximation  method  which  we  intro¬ 
duce  will  be  referred  to  as  the  Diffusion  Equation  Solution 

Sequence  or  DESS  method.  It  can  be  formulated  whenever  the 

coefficients  a(x),  b(x)  and  c(x)  are  sufficiently  differentiable 

and  the  expansions  can  be  explicitly  determined  when  the  initial 
and  boundary  data  for  the  problem  are  expressed  as  the  sum 
of  smooth  functions  and  jump  functions  (e.g.,  piecewise  poly¬ 

nomials)  . 

In  his  doctoral  dissertation  the  author1  conducted  a  rigorous 
mathematical  analysis  of  the  DESS  method  and  established 

conditions  for  the  validity  of  the  expansions  which  it  generates. 
The  problem  was  formulated  in  this  previous  work  in  a  somewhat 
different  but  equivalent  fashion  in  which  time  was  not  assumed 
to  be  small  but  a  small  parameter,  e  >  o,  multiplied  the  right 

hand  side  of  (1.1).  The  objective  of  the  current  report  is  to 
make  the  DESS  method  more  accessible  on  the  practical  level  by 
presenting  its  mechanics  in  a  straightforward  manner,  stating 
the  principal  results  concerning  its  validity  without  proof  and 
avoiding  the  mathematical  niceties  as  far  as  possible.  Our 

approach  will  thus  be  a  purely  formal  or  "engineering"  one  in 
which  (1)  the  expansions  are  sought  in  a  certain  form,  (2)  a 

necessary  system  of  equations  is  derived  for  the  individual  terms 
of  the  expansions,  and  (3)  a  method  of  constructing  explicit  solu¬ 
tions  for  these  equations  is  developed.  It  is  the  sequence  of 

functions  which  results  from  this  process  that  will  be  referred 
to  as  a  DESS.  In  the  present  work  we  shall  consider  only  the 

Cauchy  problem  and  an  initial-boundary  value  problem  with 

Dirichelet  (specified  function  values)  boundary  conditions. 
However,  the  method  is  more  general  and  can  also  be  used 
for  Neumann's  (flux  specified)  and  Robin's  (convective  heat 

transfer)  type  of  boundary  data.  The  extension  of  steps  (1) 
and  (2)  to  such  problems  is  very  straightforward  while  step 

(3)  would  require  more  effort.  In  a  forthcoming  report  we 

shall  apply  this  technique  to  the  problem  of  heat  transfer  in 
gun  barrels. 

1  J.  F.  Polk,  " Asymptotic  Expansions  for  the  Solutions  of  Parabolic 
Differential  Equations  with  a  Small  Parameter" ,  Ph.D.  Dissertation , 
Department  of  Mathematics ,  University  of  Delaware ,  Newark,  DE,  1979. 
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Let  us  now  briefly  outline  the  course  which  our  investiga¬ 
tion  will  follow.  In  Section  II  we  attempt  to  construct  formal 
asymptotic  expansions  for  the  solutions  of  (1.1)  and  obtain 
a  system  of  equations  for  the  individual  terms  of  these  expan¬ 
sions.  Then  in  the  next  three  sections  we  develop  a  number 

of  special  functions  for  the  purpose  of  solving  these  equations. 
In  Section  111  we  recapitulate  the  functions  and  H*  discussed 
in  earlier  work  and  introduce  the  new  functions  H  n  and 
H*  n  .  Next  in  Section  IV  we  define  the  functions  P-y,n  and 

Q  ’  which  satisfy  inhomogeneous  forms  of  the  diffusion  equa¬ 
tion?  Then  in  Section  V  we  formally  define  a  DESS  agd  construct 
as  particular  examples  the  sequences  (Evn>,  n  *  * 

{F  n>  and  {F  y  n)  whose  terms  are  actually  linear  combinations 
of  the  functions  Hy  _.  In  Section  VI  we  resume  the  discussion 
from  Section  II  ana  write  the  desired  asymptotic  expansions  in 

explicit  form.  Finally  in  Sections  VII  and  VI II  we  put  these 
expansions  to  use  in  constructing  approximate  solutions  for 
the  Cauchy  and  first  initial-boundary  value  problem  for  equation 
(1.1). 

This  concludes  our  introductory  remarks.  Before  proceeding 
to  the  main  discussion  we  should  first  indicate  some  of  the 

notational  conventions  which  will  be  used  in  the  following.  From 
set  theory  we  have 

R  =  all  real  numbers 

R  2  =  two-dimensional  space,  coordinates  x  and  t 
(a,b)  =  x  c  R  such  that  a<x<b 
[a,b]  =  x  e  R  such  that  a<.x£b 

H  =  R  x  (o,  •)  =  upper  half  plane,  t>  o 
Q  =  first  quadrant  in  Rl  =  (o,  •  )  x  (o,  •  ) 

Q*  =  second  quadrant  in  R2  = 

(-GO  ,o)  X  (o,  <*>  ) 

The  closure  of  a  set  in  R  or  R2  is  indicated  by  an  overbar, 
e.g. 

TT  =  R  x  [o,  ®  ) 

Q  =  [o,  ®  )  x  [o,  ® ) 

A  compact  set  S  C  R  is  a  closed  and  bounded  subset  of  R  which, 
for  practical  purposes,  may  be  considered  to  have  the  form 
[a,b]  with  a<b. 

The  derivatives  of  a  function  f(x)  will  be  denoted  by 
f  (n)(x)  and  the  partial  derivatives  of  a  function  f(x,t)  of 
two  variables  will  be  denoted 


D,f,  D  f,  D2  f 
t  x  x 

etc.  A  function  f(x)  will  be  said  to  be  of  class  C^(l)  if  it 

has  N  continuous  derivatives  in  the  open  set  I.  The  concept 
of  Hdlder  continuity  will  be  useful  in  stating  some  of  the  error 
bounds  on  our  approximations  but  will  not  be  essential  for 

their  derivation.  A  function  f(x)  is  said  to  be  Hblder  continuous 
(exponent  a,  o  <  a  <  1 )  in  an  open  set  I  if  for  each 
subset  SCI  there  exists  a  constant  K  such  that 

|f(x)  -  f(y) |  <  K|x-y|a 

for  all  x,  y  e  S.  This  property  is  denoted  by  f  e  Ca  (1); 

for  non-negative  integers  N  we  write  f  e  C  (I)  to  indicate 
that  f  £  CaU)  for  each  n  =  0,1,..., N.  Readers  not  familiar 
with  Holder  continuity  need  only  consider  the  condition 

f  e  CN+ot  (1)  as  an  intermediate  property,  stronger  than 
f  e  CN  (I)  but  weaker  than  f  e  C”1*-  (I). 

In  place  of  the  usual  Gamma  function  we  have  found  it 
convenient  to  use  the  factorial  notation 


This  is  defined  for  a>-l  by 


r  (a+1) 


oo  a  — t  . 

/  t3  e  dt 

o 


and  by  the  recursive  formula 


a!  = 


(a+1) ! 
(a+1) 


for  a<-l,  provided  a  ^  -1,  -2,  -3,...  The  factorial  function 
becomes  infinite  for  negative  integers  but  its  reciprocal  is 
well  defined  (in  fact,  entire)  for  all  choices  of  a,  and  in 
particular 


-1,  -2, 


(1.2) 


Finally,  we  wish  to  note  two  algebraic  identities  which 
will  be  used  in  Section  V: 


n  k 

L  I 

k=0  jSo 


I  I 

i=0  j  =0 


n  2  (n-i ) 

l  l  3 

j=  1  i=0  l] 


n  k  2n-l  2n-k 

k=l  j=l  k=n+l  j  =  l 


\mwm  M.IWMJIWWA'HPI 


A  proof  of  these  identities  is  given  in  the  introduction  to  refer¬ 
ence  1  and  will  not  be  included  here.  If  interested,  the  reader 
can  also  verify  these  by  writing  a  few  terms  of  the  summands 
in  a  geometric  array  and  observing  how  the  different  sides 
of  the  identities  are  obtained  by  summing  the  array  elements 
in  different  orders. 


II.  FORMAL  EXPANSION  PROCEDURES 

The  parabolic  differential  equation 

Dtu  =  Lu  s  a(x)D^u  +  b(x)Dxu  +  c(x)u  (2.1) 

is  of  fundamental  importance  to  transport  processes.  The  un¬ 
known  function  u(x,t)  can  be  thought  of  as  any  transport  quanti¬ 
ty,  such  as  temperature,  species  concentration,  vorticity,  etc. 
and  the  terms  on  the  right  hand  side  of  (2.1)  correspond  to 
contributions  to  its  rate  of  change  due  to  diffusion,  convection 
(advection)  and  species  production  (or  consumption)  respectively. 
In  the  present  form  the  coefficients  depend  only  on  the  spatial 
coordinate  x  and  thus  equation  (2.1)  is  linear  although  it 
is  frequently  non-linear  in  actuality.  Nevertheless  this  equation 
is  pertinent  to  non-linear  problems  for  two  reasons:  (1)  a 

proper  understanding  of  non-linear  processes  can  only  be  achiev¬ 
ed  after  a  thorough  rendering  of  the  linear  case,  and  (2)  a 
non-linear  equation  can  be  solved  numerically  by  quasi-lineariza¬ 
tion  in  which  non-linear  coefficients  are  expressed  as  functions 
of  the  spatial  variable  only  during  a  single  time  step.  This 
procedure  is  valid  for  sufficiently  short  time  intervals  and 
can  be  performed  repeatedly  to  advance  the  computation  over 
longer  time  spans.  In  computational  fluid  dynamics  even  simpler 
versions  of  (2.1)  are  sometimes  studied2  to  gain  insight  into 
the  nature  of  transport  processes.  In  fact,  it  seems  that  the 

majority  of  the  working  tools  used  in  practice  (stability  criteria, 
order  of  magnitude  error  estimates,  etc.)  have  been  rigorous¬ 
ly  justified  only  in  the  context  of  linear  initial  value  problems. 

In  order  to  motivate  more  complicated  expansion  procedures 
for  obtaining  approximate  solutions  of  (2.1)  let  us  first  intro¬ 

duce  the  simplest,  most  direct  method.  We  suppose  that  initial 
data  has  been  specified  in  the  form 

u(x,0)  =  f(x)  (2.2) 

valid  for  x  in  a  given  interval  I  C  R.  The  N-term  Taylor 

expansion  for  u  in  the  time  direction  from  the  initial  line  is 

then 

2  P.  J.  Roaahe ,  ”1 Computational  Fluid  Dynamics",  2nd  Edition,  Hermosa 
Publishers,  Albuquerque,  NM,  1976. 


r 


* 


N  (k )  u 

l  DjK;  u ( x ,0  )  tK/k! 

k=0 

But,  by  repeated 
as 


or 

N  /  u  \  u 

l  L'K  f ( x )  t/k ! 

k=0 
(k) 

where  L  denotes  the  operator  L  applied  k  times.  Because 

of  its  special  significance  we  shall  refer  to  the  last  expression 
as  the  N-term  regular  expansion  for  u  and  denote  it  by 


formal  substitution  of  (2.1)  this  can  be  written 

N  /i.  \  k 

l  L'K'  u(x,0)  tK/k! 
k=0 


U 


N 


N 


L(k)  f(x)  tk/k ! 


(2.3) 
(k) 


(x,t)  =  J 
k  =0 

This  expansion  is  well  defined  if  all  of  the  terms  L'1''  f(x) 
are  meaningful;  this  holds  when  a,  b,  and  c  have  2  N  ~2  deriva¬ 
tives  and  f  has  2N  derivatives  at  x.  Using  slightly  stronger 
assumptions  than  these,  the  regular  expansion  has  been  shown 
to  accurately  approximate  u  and  we  shall  state  these  results 
formally  in  Sections  Vll  and  VI 1 1 . 


Assuming  that  a,  b  and  c  are  suitably  behaved,  the  regular 
expansion  breaks  down  in  two  not  infrequent  situations: 

(1)  The  initial  value  function  f(x)  is  not  sufficiently  smooth 
and  thus  U  is  not  even  defined.  Since  f(x)  is  usually  determin¬ 
ed  by  physical  measurement  or  inference,  there  is  no  guarantee 
that  it  will  have  any  particular  "analytic"  form.  For  example, 
bringing  both  hot  and  cold  objects  suddenly  together  leads 
to  discontinuous  initial  values.  On  the  other  hand  it  can 
be  supposed  that  almost  any  physically  reasonable  f(x)  is 

at  least  piecewise  smooth  with  jumps  of  various  orders  occuring 
at  discrete  points.  For  this  reason  a  particularly  relevant 

choice  of  initial  values  is  the  jump  function 

f(x)  =  h y  (x-xq)  (2.4) 

where 

hY (z)  =  {  zVy!  z  >  0  (2.5) 

0  z  <  0 


with  y>0. 
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(2)  a  boundary  condition  can  be  imposed  along  a  boundary 
such  as  x=xQ  which  is  incompatible  with  the  regular  expansion 
even  when  well-defined.  In  fact,  the  most  that  can  be  expected 
in  general  is  that  the  boundary  condition  is  continuous  with 
the  initial  data  at  the  corner  point  (xQ,o).  This  would  ensure 
that  only  the  lowest  order,  but  no  higher  term  of  the  regular 
expansion  agrees  with  the  prescribed  boundary  data.  Such 
incompatibilities  are  evidenced  by  "thermal  boundary  layers" 
in  short  term  heat  conduction  problems.  The  most  fundamental 
boundary  condition  is  that  in  which  the  function  values  are 
specified  (Dirichelet  condition)  along  a  fixed  boundary  and, 
among  these,  the  most  important  are  those  which  are  homogeneous 
or  can  be  expressed  as  a  power  of  t. 

From  these  considerations  we  are  led  to  formulate, 
as  "canonical"  problems,  equation  (2.1)  together  with  the  follow¬ 
ing  choices  of  initial/boundary  data 

(IVP)  u(x,o)  =  hv  (x-x  )  xeR 

o  Y  o 

( BVP )  u(x,o)  =0  x  <x 

o  =  o 

u(xQ,t)  =  h,^  o<t,<At 

Where  y>0  and  At>0.  (The  use  of  y/2  rather  than  y  in  (BVP)0 
leads  to  more  conveniently  stated  results.)  The  upper  bound 


on  time 

is  indicated  by 

At  to  emphasize  our 

interest  in  short 

duration 

problems 

approximations. 

Let  us  also  formulate 

the  two  converse 

(1VP)J 

u(x,o) 

=  h*  (x-xQ) 

X  e  R 

where 

h*(z) 

=  h^  (-z),  and 

(BVP)# 

u(x,o) 

=  0 

x±xo 

o 

u(xo,t) 

=  h r/2  (t) 

0<t<  At 

The  expansions  for  these  latter  two  problems  can  be  obtained 
from  the  preceding  ones  by  a  straightforward  transformation. 
We  shall  state  the  parallel  results  for  convenience  without  addi¬ 
tional  comment. 

From  the  theory  of  parabolic  differential  equations’  4  5  it  is 

3  A.  Friedman ,  "Partial  Differential  Equations  of  Parabolic  Type ", 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ,  1964. 

4  0.  A.  Ladyzenskaja,  V.  A.  Salonnikov,  N.  N.  Ural'ceva,  "Linear 
and  Quasilinear  Equations  of  Parabolic  Type",  American  Mathematical 
Society  Translations  of  Mathematical  Monographs,  Volume  23,  1968. 

5  P.  C.  Rosenbloom,  "Linear  Partial  Differential  Equations",  Surveys 
in  Applied  Mathematics  V,  John  Wiley  and  Sons,  Inc.,  New  York,  1958. 
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well  known  that  a  solution  u(x,t)  exists  for  each  of  the  prob¬ 
lems  just  formulated;  furthermore  it  will  be  unique  if  restricted 
to  the  class  of  functions  such  that 

|u(x,t)  |  <  K  exp  (k  x2)  (2.6) 

for  some  constants  K,  k^O,  uniformly  for  x,t  in  the  domain  of 
interest. 

Let  us  now  proceed  with  the  formal  development  of  asymptotic 
expansions  for  u(x,t)  basing  our  derivation  on  two  considerations: 
first  our  interest  is  primarily  in  the  local  region  near  x=x 
(since  the  regular  expansion  suffices  elsewhere)  and  second0 
we  are  concerned  only  with  the  short  time  behavior.  We  shall 
suppose  that  equation  (2.1)  has  already  been  written  in  a  form 
such  that  u,  x  and  t  are  non-dimensional  with  t  restricted  to 
the  interval  o<t£At  where  At<<l. 

To  emphasize  the  local  short  term  behavior  of  u  let  us  intro¬ 
duce  the  following  "stretched"  variables 


o  =  (x-xQ)/e 
t  =  t/e2 


(2.7) 


where,  for  convenience,  we  henceforth  let  e  =  /At. 

In  terms  of  these  variables  we  may  re-express  the  solution  u 
as  the  transformed  function  u  defined  by 


u(a,T)  =  u(xQ+  eo,  e  2 x )  =  u ( x , t ) 

Similarly  the  coefficients  a,  b,  and  c  may  be  expressed  as 

a(c )  =  a(x  +  eo) 

6 (o  )  =  b(xQ+  eo) 

C(  o)  =  c(x  Q+  eo) 

Substitution  of  these  into  (2.1)  leads  to  the  following  equation 
for  u 


D  u  =  a(o)D2u  +  e  6(o)D  u  +  e2  c(o)u  (2.8) 

t  o  o 

In  this  form  standard  perturbation  methods  can  be  employed 
to  develop  expansions  for  u  in  terms  of  the  small  parameter  /aT. 
Let  us  suppose  that  u  can  be  expanded  in  the  form 

U  -  ep  [u  +  e  U,+  e2U,+  ...  (2.9) 

0  1 

where  p  is  not  yet  specified  but  depends  on  the  particular  choice 
of  boundary/initial  data.  If  the  coefficients  a,  b  and  c  are 
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sufficiently  differentiable 
term  Taylor  expansions, 
o  these  are 

at  x=x0  then 
In  terms  of 

we  can  form 
the  stretched 

their  n- 
variable 

a(  o)  = 

l  a,  (e  °)  /k! 
k=0  K 

+  A 

n 

6(a)  = 

f  b.  (c  )  k/k ! 

k=0  K 

4  B 

n 

c{  a )  = 

f  c.  (e  o)k/k! 
k=0  K 

4-  C 

n 

where  a^,  b^ 

and  c,  denote  the  Taylor  coefficients 

a  =  a(k)(x  ) 
k  o 

,  ,(k),  v 

b,  =  b  ( x  ) 

k  o 

(k) ,  . 

ck  =  c  (xo) 

(2.10) 

and  A  ,  Bn  and  Cn  are  remainder  terms.  If  (2.9)  -  (2.10) 
are  substituted  into  (2.8)  and  terms  with  similar  powers  of 
e  collected  together  then  the  following  system  of  equations  results 


and  in  general 


-  aoD'„)  ao  ■  0 
ID  -  aJ>’J  5i  * 


O  O 


I  u0 


[Dt  "  a0D03  G2  =  ^iGi  +  L2  G 


[D  T  -  a  D 2  ]  u  =  7  L.  u  , 

L  t  o  oJ  n  .  “ .  k  n-k 
k=l 

for  n=l,  2,  ...  where 

k  k-1  _  0k-2 

-  u. 


L,  u 
k 


a,—  D2u  +  b,  - 

kk  !  k_1  (k-1 ) ! 


D„u  +  c,  0 
a  k-2 


(k-2) ! 


(2.11) 


(In  this  notation  we  are  making  use  of  (1.2)  in  the  case  k=l.).  • 
The  operator  is  clearly  well  defined  whenever  a,  b,  and 

c  have  k,  k-1  and  k-2  derivatives  at  x0,  respectively. 

Next,  by  rewriting  the  boundary  and  initial  data  for  ( 1  VP ) 0 
and  (BVP)0  in  terms  of  the  stretched  variables  we  obtain 

u( o,o)  =  h  (o  )  o  e  R 

for  problem  (IVP)  ,  and 

u(o ,o)  =  0 


o  <  o 


\  W  ,  l  / 


(t) 


for  problem  (BVP)0.  From  this  it  follows  that 

.  y  in  problem  ( 1  VP ) 0 
”  -  1  0  in  problem  (BVP)q 


o<t  <  1 


(2.12) 


and  that  the  individual  terms  of  the  expansion  for  u  must  satisfy 
the  supplementary  conditions 


Gk(a’o)={0 
for  problem  (1VP)  and 

°  uk(o,0)  =  0 

u,  (o, t)  =  J  ^ 
0 

with  o,  t:>0  for  problem  (BVP)  . 


,(° 


y/2 


)  for  k=0 
for  k=l,  2,  ... 

for  k-0,1,2... 

(  t  )  for  k=0 

for  k=l,2,... 


(2.13) 

(2.14) 

(2.15) 


The  forgoing  systems  of  equations  and  initial/boundary 
conditions  for  the  terms  u,  have  a  unique  solution  among  the 
class  of  functions  satisfying  (2.6)  which  can  be  uniquely  deter¬ 
mined.  To  accomplish  this  we  suspend  our  present  discussion 
for  the  moment  and  in  the  following  three  sections  investigate 
some  relevant  special  functions.  Readers  not  interested  in  the 
details  can  continue  in  Section  VI. 


III.  THE  FUNCTIONS  H  AND  H 

Y  Y ,  n 

The  functions  have  been  discussed  in  previous  work1'  6  . 
For  convenience  we  briefly  recall  their  definitions  and  basic 
properties.  Wheny>-1  we  define 

00 

H  (x,t)  =  (4nt)  ^  f  (sy/y!)  exp  (-(x-s) 2  /4t)ds  (3.1) 

1  *o 

for  t  >0 ;  for  y£-l  H  Y  is  defined  recursively  by 

HY(x,t)  =  Dx  HY+1(x,t)  (3.2) 

for  t  >0.  (This  also  holds  as  an  identity  for  y  >  —1 . )  Along 
x=0  these  functions  are  assigned  the  values 

H  y(  x,o)  =  h^  (x)  xeR  (3.3) 

where  hy(x)  is  defined  by  (2.5);  except  at  the  origin  x=t=0 
with  y  <6  this  can  be  seen1  to  form  a  continuous  extension  of  (3.1) 
and  (3.2).  These  functions  satisfy  the  relationships 

6  J.  F.  Folk,  "Special  Function  Solutions  of  the  Diffusion  Equation", 
ARBRL-TR~02182,  US  Army  Ballistic  Research  Laboratory,  1979. 

(AD  #A075324) 
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YHy(x,t)  =  x  Hy_1  (x,t)  +  2t  Hy_2  (x,t) 


Hy(x,t)  = 


Hv  (x  a  t) 


(3.4) 

(3.5) 


for  any  y  eR  and  a>0.  Particular  examples  of  Hy  functions 
are 


where 


H0(x,t)  =  (1/2)  erfc  (-x/ 

H_j(x,t)  =  Uirt)-*^  exp  (~x2/4t) 
H_2 ( x , t )  =  ~(x/2t)  H_j  (x,t) 

erfc(z)  =  (4/ir)^  /  “exp  (-s2 )ds 

» 


(3.6) 


The  reflected  functions  h*  and  H*  are  defined  by 

Y  Y 


h*  (z)  =  hy(-z) 

Hy  (x,t)  =  Hy  (-x,t) 


(3.7) 

(3.8) 


All  of  the  functions  Hy  and  H*  were  seen  to  be  infinitely  differ¬ 
entiable  with  respect  to  x  and  t  for  t>0  and  to  satisfy  the 


formulas 


satisfy 


Dx  Hy  "  Hy-1 
Px  H*y  -  -H*~1 


Dt  Hy  -  Hy -2 
Dt  Hy  "  Hy-2 


9 


Thus  they  are  solutions  of  the  heat  or  diffusion  equation 

[D,  -  D*]  u  =  0 
t  x 

for  t>0.  Along  x=0  they  assume  the  values 


Let  us  now  define 


Hy'o,  t)  =  H*  (o,t)  =  /Ty/2(y/2)! 


Hy,n  (x’°  =  XT  Hy  {x’t} 


H*  (x,t)  =  H  ( —x  ft)  = 
Y  *n  y » n 


“I —  H  (-x,t) 
n !  y 


(3.9) 


(3.9)’ 


(3.10) 


(3.11) 


(3.12) 


(3.12)* 


where  y  e  R  and  n  is  an  integer.  These  can  be  regarded  as 
generalizations  of  the  functions  H  y  since  setting  n=0  yields 


H  Y,o(x>t)  *  Hy(x.t) 


(3.13) 


(3.13)* 


H*0(x,t)  *  H*(x,t) 

Moreover,  for  negative  integers  n,  these  functions  are  seen  to 
vanish  identically 


H 


Y.n 


H*  =  0 

Y.n 


(3.14) 


ri=-l,-2,-3,  ...  ,  in  view  of  (1.2).  Using  (3.3)  and  (3.11) 

we  see  that  along  t=0  and  x=0  these  functions  take  on  the  values 


H 


Y,  n 


(x,o) 


H  (o.t) 

Y.n 


=  hy(x)  hn(x) 

{’Vt) 


if 


Y+n 


n=0 


T 


x>0 
x<  0 


if  otherwise 


(3.15) 

(3.16) 


From  differentiation  formulae  (3.9)  and  Leibnitz's  rule  for  repeat¬ 
ed  differentiation  of  a  product  we  obtain 


and 


D  H 
x  Y.n 

Dk  H 
x  Y,n 

D.  H 

t  Y.n 


HY,n-l  +  HY-l,n 

r  (k)  H  . 

j4_0  J  Y-J.n-k+] 

H  . 

Y-2,n 


[D2  -  Dj  H  „  =  2  H  ,  ,  +  H  n  „ 

x  t  Y«n  y  — 1 » n— 1  Y.n-2 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


For  any  constant  a>0  it  then  follows  that  Hy  n  (x,at)  satisfies 
the  following  inhomogeneous  form  of  the  diffusion  equation 


[aD2-D  ]  H  (x.at)  =  2aH  ,  .(x,at)  +  aH  0(x,at)  (3-21) 

for  t>0.  Also,  for  any  a>0,  we  can  obtain  the  following  identity 
from  (3.5) 

Hyn(x,t)  =  a'(Y+n)/2  H  y  (x/a,at)  (3.22) 

In  Reference  1  these  functions  were  studied  in  greater  detail, 
in  particular,  their  asymptotic  behavior  for  large  |  x|  was  char¬ 
acterized.  In  the  present  work  we  do  not  require  this  analysis 
except  to  recall  that 

|H  (x,t)  |  <_  const,  [h  (x)  +  /t  Y+n]  (3.23) 

Y*  n  Y  +n 

and  thus  H  n  clearly  belongs  to  the  class  of  functions  satis¬ 
fying  (2.6).  ’ 
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IV.  THE  FUNCTIONS  P 


These  functions  are  thus  linear  combinations  of  the  n  func¬ 
tions.  Their  significance  is  indicated  in  the  following  ’two  pro¬ 
positions.  Recall  out*  notation  Q  and  Q*  for  the  first  and  second 
quadrants  in  R  2. 

If  y+n>-2  then  (4.1)  defines  the  unique  function 
which  is  continuous  in  §*,  uniformly  bounded  whenever  t  is 
bounded,  which  satisfies  the  equation 

1°.  -  Dl*  <3y,n(l‘',)  '  **Y  ,n*x,t) 

in  Q*  and  takes  on  the  homogeneous  boundary  and  initial  values 


c 

*> 

>- 

O’ 

(x,o) 

=  0 

x^O 

(4.5) 

Q 

x  Y.n 

(o,t) 

=  0 

t>0 

(4.6) 

iasuudl  The  uniqueness 

of 

such 

functions  is  well 

known  and 

follows,  for  example,  from  the  results  of  Widder.7  The  continuity 
of  0  Y,n  for  t>0  is  established  in  Proposition  (1.3.1)  of  Reference 
1.  (The  restriction  Y+n>-2  is  required  here.)  The  homogeneous 
initial  and  boundary  values  follow  directly  from  substitution 
of  (3.15)  and  (3.16)  into  (4.1).  The  boundedness  requirement 

7  D.  V.  Widder ,  " The  Heat  Equation ",  p-139}  Academic  Press ,  New  York , 
1975. 


is  directly  verified  by  applying  inequality  (3.23)  to  the 
individual  terms  of  Q  .  Thus  it  only  remains  to  show  that 

Q  v  _  satisfies  (4.4).  From  (3.20)  we  have 

T » n 


n+1 


[D«  -  Di1Qr.n  *  -J/-5 


1  ,n+2-k  r ,  u  1 

Y+n+l-k,k-l  Y+n+2-k,  k-2 


n+1 

L 


,  hn+l-k  „  r.  ,  l,n+2-k  H 

1  Y+n+l-k,  k-1  ^2  "2  Y+n+2-k,k-2 


where  the  k=l  term  in  the  second  sum  was  omitted  in  view  of 
(3.14).  Shifting  the  index  in  this  sum  by  setting  j  =  k-1  we 
have 

n+1  .  ,  . 

lt,t  -  DxlQ„n  ‘  l-7>  H 

-  ?  (-  -i)n+1-J  H 

jSl  1  1 

=  H 


Y+n+l-k,  k-1 

^Y+n+l-j,  j-1 


Y*n 

This  completes  the  proof. 


If  Y +n  >  -2  then  (4.2)  defines  the  unique  function 
which  is  continuous  for  t>0,  satisfies  the  equation 


iDt  '  D0  *  N.n1*'" 

in  H,  with  initial  values 


(4.7) 


P  (x,o)  =  0 

Y  ,n 


xE  R 


(4.8) 


and  which  is  bounded  such  that  for  any  T>  0  constants  k, 
K>0  exist  for  which 


uniformly  for  o<t<T 


P  (x,t)|<_K  exp  (kx2  ) 
'  *  n 


PROOF: 


Essentially  the  same  as  for  Proposition  (4.1)  except 
that  (4.8)  must  be  established  for  x>0.  But  the  constant  Cy>n 
was  chosen  precisely  so  that  this  result  holds,  as  may  be  seen 
by  substituting  (3.15)  into  (4.2)  with  t=0. 


_ _  If  we  denote  by  H 

linear  combinations  of  functions  H 


RFMAI-'K 


x  the  family  of  all  finite 
with  y  +n=  x  then  clearly 


Y,n 


wiiiu  111a,  uwno  ui  iuulliwuo  1 1  y  .  wun  I  th-  a  men  v-icouy 

and  QY>n  belong  to  HY+n  .  ’The  importance  of  the  above 
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propositions  is  that  for  any  v  in  the  solution  to 

[D.  -  D* ]  u  =  v 
t  x 

with  homogeneous  initial/boundary  data  can  be  found  in  the  class 

Hx+2*  It  is  this  behavior  which  enables  us  to  obtain  solutions 

for  the  equations  of  Section  11  explicitly  in  terms  of  the  Ry  n 

functions.  In  fact,  for  the  canonical  problems  (IVP)Q  and  (B^P)0, 

the  term  of  the  asymptotic  expansion  (2.9)  will  be  found 

in  the  class  H  . 

r+k 

For  any  a>0  the  foregoing  results  imply  that 

the  functions  Pv  (x,at)  and  Qv  „  (x,at)  solve  the  equation 
T  i  n  t  t  n 

[D,.  -  aD2  ]  u(x,t)  =  aH  (x,at) 
t  x 

in  place  of  (4.4)  and  (4.7)  and  satisfy  the  same  homogeneous 
initial/boundary  data  as  before. 


IRHMARK 


The  functions 


and 

Q 

satisfy  the  equations 


p* 

Y.n 

(x,t) 

■  P,.n 

( —x ,  t ) 

Q* 

VY.n 

s 

(x,t) 

QY»n 

( —x ,  t  ) 

[Dt  - 

D2] 

X 

P*  „  = 

Y.n 

H* 

Y  . 

[Dt  - 

D2] 

X 

Q* 

Y.n 

H*  n 

Y.n 

in  H,  and 

[D,  -  D2]  Q*  =  H* 

t  xJ  ^Y,n  Y.n 

in  Q,  and  take  on  the  initial  and  boundary  values 


* 

(x.o)  =  0 

X  e  F 

Ytn 

i 

(x.o)  =  0 

x  >0 

Y,n 

(o,t)  =  0 

t>0. 

Y,n 

* * 

V.  DIFFUSION  EQUATION  SOLUTION  SEQUENCES: 

THE  FUNCTIONS  Ev  n  AND  Fv  n 

T  9  11  If11 

We  have  now  reached  the  point  at  which  we  can  formally  con¬ 
struct  solutions  for  the  equations  derived  in  Section  II. 

Let  N>0  be  an  integer  and  let  be  an  (ordered)  parameter 
set  given  by 

it n  =  (a0,ai,  ...  ,  a^,  bQ.b^,  ...  , 
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bN-r  co’  cr 


CN-2}  ' 


These  parameters  may  of  course  be  considered  as  the  Taylor 
coefficients  of  a(x),  b(x)  and  c(x)  at  x=x  but  their  genesis 
is  unimportant  here  and  the  discussion  may  pPoceed  independently 
of  subsequent  applications. 


For  n=0,l,  N  we  may  define  operators  Ln  by 


L  u  =  a  — -r  D2u+b 


,  7- — r-TT  D  u+c  „ 
n-1  (n-1 ) !  x  n-2 


u.  (5.1) 


In  particular 


L  u  =  a  D2  u 
o  ox 

L.u  =  a.xD2u+b  D  u. 

1  1  x  o  x 

A  sequence  of  functions  U  =  {u  (x,t):  n  =  0,1,.,..,  N) 
defined  in  the  closure  of  a  domain  D  C  H  will  be  called  a 
Diffusion  Equation  Solution  Sequence  (DESS)  in  D  with  respect 
to  the  parameter  set  if  the  following  equations  are  satisfied 


[Dt  - 


f° 

L  ]u  =  / 
o  n  ^n 

1  y  l,  u 


k  n-k 


if  n  =  0 


if  n  =  1,2,  ....  N 


(5.2) 


If  D  is  unbounded  we  also  require  that  for  any  T>0  there  exist 
constants  k,  K  >  0  such  that 


un(x,t)  I  1  K  exP  y2  ) 


(5.3) 


Uniformly  in  0<_t£T  for  each  n  =  0,1,  ....  N.  We  shall  use 

the  notation  Ue  P(D,  Iljj)  or,  more  simply,  UeP^  to  indicate  that 
a  sequence  U  =  (un:  n  =  0,1,  ...,N)  is  of  the  type  described. 

If  U,V  e  Vfj  and  a,  B  e  R  then  we  may  define  the  formal  sum 
sequence 

W  =  aU  +BV 

termwise  by 

w„  =  «u„  +  Bv, 
n  n  n 

( k ) 

and,  for  0  ^k  <_N,  we  may  define  the  "shifted"  sequence  U 
termwise  as  the  sequence  whose  n-th  term  is 

,  0  0<n<k-l 

Un  =  { 

n  1  in 

%_k  k!"<N. 

(The  superscript  k  here  is  not  meant  to  imply  differentiation.) 
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^5 


It  should  be  clear  that  these  sequences  are  also  of  class  V ... 

Moreover,  let  c  R  and  for  k  =  0,1,  ...,N,  with  ^1)^  = 

{U.  },  and  define  the  sequence  of  "shifted"  partial  sums  by 

k  f  n 

N  (k) 

V  =  I  d  u 

k=0  K  K 

or  termwise  by 

v  a  ,  (k) 

Vn  '  Jo  k“k’n-k  ’  Lo  kUk’n 

Then  the  previous  observation  implies  that  Vep^  . 


From  the  existence  and  uniqueness  properties  of  the  diffusion 
equation,  D^u  =  Lqu,  we  see  that  we  may  specify  the  individual 
terms  of  a  DESS°  by  requiring  that  they  also  satisfy  certain 
boundary  and  initial  data.  The  following  sequences  will  be 
of  particular  interest: 


(1)  E  and  E# 
Y.n 


Y.n 


y  ,  n 

h  (x) 

if 

n 

=  0 

(x,0)  = 

{  Y 

1  0 

if 

n 

=  1,2 

h*(x) 

if 

n 

»  0 

(x.O)  = 

{  Y 

1  0 

if 

n 

=  1,2, 

(2)  F  and  F* 

and  satisfying’  ’ 


Y.n 
'#  n 


(5.4) 

.  ,N. 

(5.4) # 

.,  N. 

defined  in  (J*  and  ^  respectively 


Y.n 


(x.O) 

=  0 

for 

xIO, 

n  =  0,1,.. 

,.,N 

(5.5) 

(x,0) 

=  0 

for 

x>0, 

n  =0,1,., 

...N 

(5.5)# 

/ry 

=  0 

(0,t) 

=  F# 

_(0» 

uJ 

T7727T  n 

(5.6) 

y 

,  n 

1 

0  n 

=  1 ,2, ...» N 

for  t>0. 


The  following  two  propositions  are  rather  technical  but 
establish  that  each  of  the  terms  of  the  foregoing  DESS’s 
is  expressible  as  a  linear  combination  of  functions  H  y  n  and 
H%  n  with  coefficients  that  can  be  explicitly  determined  from 
recursion  formulae.  The  boundedness  condition  (5.3)  will  then 
be  automatically  satisfied  due  to  the  growth  properties  of  such 
functions  [see  (3.23)]. 
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(5.7) 


Let  0  <  n  <_  N,  then 
2n 

E  (x,t)  =  T  a  ,  H  ,  ,  (x,a  t) 

Y,n  k=0  Y,n,k  Y+n_k’  k  ° 


where  the  coefficients  a 
recursive  formulae: 


Y,n,k 


for  n  =  0 


for  n  =  1,2, . . . , N  a 


Y»n,k 


where,  from  (4.3) 


'y+n-2-j,  j 


are  determined  by  the  following 

“  ^,0,0  =  1 
2n-l 

l0  CY+n-2-  j  ,]aY»n,  J 

,  .  lf  k  =.°  (5.8) 

l  (4)  “r.n.i‘k-2 

j  =  l 

if  k  =  1 ,2, . . . ,2n 


,j+2-i 

i.1!  4  m 


and  the  tr  .  are  defined  for  i  =  0,l,...,2n-l  by 
Y,n,] 

1 


*  J 

=  th  ai  ^  Y.n-i.j-i 

*  j  l  2ai  (i)  +bi-l  G-l)^ay  ,n-i,  j-i+1 


(5.9) 


*  .1  [ai(i)+bi-l(i-l)+Ci-2(i-2)]aY,n-i,j+2-i 
i=l 

with  the  upper  summation  limits  being  given  by 

1  "  { 

^n-j  li  n+i  ^  j  <-  i— x 

j+1  if  0  <  j  <  r-1 

J  =  {  (5.10) 

K  ={ 
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j 

if 

0 

1  j 

< 

n 

2n-j 

if 

n+1  < 

j 

<  2n-l 

j+1 

if 

0 

1  j 

< 

n-1 

2n-j-l 

if 

n 

1  j 

< 

2n-l 

j+2 

if 

0 

<_j 

< 

n-2 

CM 

1 

1 

C 

CM 

if 

n- 

1  < 

i  . 

<  2n-l 

It  should  be  understood  that  the  sums  appearing  in  (5.9)  are 
omitted  when  the  lower  summation  limit  exceeds  the  upper. 


The  proposition  holds  for  n  =  0  since  we  then  have 
E  (x.t)  =  H  (x,a  t) 

Y.O  Yi°  O 

Proceeding  by  induction  we  may  assume  that  (5.7)  is  valid 
for  n  =  0,1,..., m-1  where  1  <  n  <_  N  and  show  that  is  must 
also  hold  for  n  =  m.  For  convenience  in  this  proof  we  shall 
write  Hy  n  for  HY>n(x,a  t).  Hie  k-th  term  on  the  right  hand 
side  of  (5.2),  when  un  =  E  and  n  =  m  is  then 


k  2{m-k) 


L.  E  ,  =  a, 

k  Y,m-k  k 


a.  rr  T  a  .  .DSH  .  .  . 

k  k!  Y.m-k,]  x  Y+m_k-J,J 

^k-1  2(m-k) 

+  bk-l  (k-l)T  “y  ,m-k,jDxHY+m-k-j,  j 


k-2  2(m-k) 


+  ck-2  (k-2)T  aY,m-k,  jHY+m-k-j,j 

where  1  <_  k  <_  m.  Using  (3.17)  and  (3.18)  yields 

k  2(m-k) 

bk  Y»m-k  k  Z0  °  Y»nt-k,  Y+m-2-k-j,j 


+2Hy+m-l-k-j ,  j-1  +HY-Hifk- j ,  j-2  V 
k-1  2(m-k) 


3k-l  (k-1 ) !  aY,m-k,  j^HY+m-l-k-j,  j+HY +m-k-j,j-l  ^ 


k-2  2(m-k) 


c-2  Ck-2)!  .  I, 


aY,m-k,j  Y+m-k-j,j* 


Using  the  identity 


H  (x.t)  =  (m+n)  H  (x.t) 

m !  y  * n  n  y  » m+n 


which  follows  directly  from  definition  (3.12)  we  then  have 
2(  m-k ) 

bk^Ytin-k  "  .  I  ak  Ik  J  a y, m-k, i^Y+m-2-j-k, j+k 

]  =o 
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2(m-k)  .  .  , 

*  IQ  (2.k  (Jf) 


^k-l^la  y , m-k , j 


Y+m-l-j-k,]+k-l 

2 (m-k)  ,  _ 

*  jlo  [a“  (  ^  +bk-l 


fj+k-2)  +C  ^+k~21l 

lk-1  J  ck-2  U-2 


•  H 

'aY,m-k,j  Y+m-j-k, j+k-2 

The  right  hand  side  of  (5.2)  when  u  =  EY  and  n  =  m  is 
obtained  by  summing  this  expression  with  respe’ct  to  k  =  l,2,...,m. 
For  m  =  1  this  reduces  to 

L.Ev  „  =  a.H  0  1+b  H  . 

1  Y,o  1  y-2,1  O  Y-l.O 

For  m  ^  2  we  may  reexpress  the  summation  with  the  help  of 

identity  (1.4);  thus 


m  m  n 

k  =  l  ai  ^  aY.ni-i,n-iHY+m-2-n,n 


2m- 1  2m-n 

+  ^  Si  (  i )  aY.m-i,n-iHY+m-2-n,n 


n=m+l 
m  n 


+  T  Y  [2a.  f11:1!  +b.  .  (  n*!)]a 

n  -  1  i=l  1  1  1  1  '  i-l  1  i— 1'**  y  * m_l » n~l 


n-K 


•  H  , 

Y+m-l-n,n-l 

2m-l  2m-n  ,  . 

+  y  i  i 2a.  C1:1)  +b.  .  .  . 

nim+1  1  =  1  11  1-1  x~l  Y’m-i.n-i 


•  H 


Y+m-l-n,n-l 
m  n 


*  nil  L  ("‘2)  *b‘->  ^ 


■a  H 

Y^m-i.n-t  Y+m-n»n-2 
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2m- 1  2m-n  „  „  „ 

♦  i .  .1.  •. inflow  (?:?>  «i-2(d)i 

n=m+l  i=l 


•  ft  _  1 1  n  . 

Y,m-i,  n-i  Y+m-n»n-2 

Shifting  indices  and  rearranging  slightly  we  get 

m  n 

,m-k  ,1.  .1  ai  (  i  )  ay,m-i,n-i^ 


(5.11) 


>  L.  E 
kil  k  Y 


m-1  n+1 


*  *  * '  t » •••  ****  *  Y +m  2  n,n 

2m- 1  2m-n  „  (5.12) 

+  T  y  a.  fn)a  .  .H  „ 

n-m+1  i-1  1  Y* m_1 » n_l  Y+m-2-n,n 


+  n=0  i-1  [2&i  ^  +bi_1  ^  i-l)*  aY.m-i,n+l-i  HY+m-2-n,n 

2m-2  2m-n-l 

+  T  y  [2a.  (n)  +b.  .  f  n,)]<x  1  •  H  „ 

^  1  l  lu  l-l  ti-yJ  Y.m-i,  n+l-i  Y+m-2-n,n 

m-2  n+2 

+  y  y  ra.  fni+b.  .  r.H.)  +c.  0  f  n„iia  .  „  .  h 

n_0  i-i  1  1  i-l  u-ll  i-2  1  i-2<J  y» m— i > n+2— i  Y+m-2-n,n 

2m-3  2m-n-2 

+  E  .  1.  tai(i)  +bi-i  C i— 1  )+ci-2  (i-2^®  Y,m-i,n+2-iHY+m-2-n,n 

n=m-i  i*i  11 

Collecting  together  the  coefficients  of  similar  terms  we  can 
then  write 


m 

y  L..E 


2m -1 


k^=l  *'k‘'Y,m-k  ”  J=0  a  y  .m,nbY+m-2-n,n 


(5.13) 


where  the  a  are  given  by  (5*9)  above  for  m  _>  1.  In 

the  nt=l  Y,n^rPost  of  the  terms  are  absent  leaving  only 

oT  ,  n  =  b  and  ST  ,  ,  =  a. . 

Y.l.U  o  Y  *1.1  1 

Equation  (5.2)  applied  to  the  m-th  term  of  the  DESS  {EYn) 
can  thus  be  written  as 

2m-l 

[D.  -  a  D*]E  (x,t)  =  7  g  H  0  (x,a  t). 

t  o  x  y  *  m  ^_q  y  » m  * n  y  +m-2-n ,  n  o 

Recalling  Proposition  (4.2)  and  Remark  (4.2)  and  using  linear¬ 
ity  it  immediately  follows  that  for  m  =  1,2... 
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ft 


E  m(x,t) 
y»  in 


_1_ 

a 


a 


2m-l 

7  "a  P  0  (x,a  t) 

n=0  Y*m*n  Y+m-2-n,n  o 


2m-l 


n+1 


E  a  [  7  f— i)  H  ,  ,  ( x ,  a  t ) 

n=0  Y,m,n  k=l  y+m-k.k  o 


n+2-k 


-  c  „  H  (x,a  t)  l 
Y+m-2-n,n  Y+m  o  1 

By  collecting  coefficients  with  the  help  of  (1.3)  this  can  be 
rewritten  as 

2m 

E  (x,t)  =  7  a  i  H  ,  .  (x,a  t) 

Y.m  k^0  Y*m»k  Y-fm*k»k  ° 

where  “„fm,k  is  given  by  (5.8)  with  n  =  m  M.  This  completes 
the  pi  oof. 

For  n  =  0,1,2  the  coefficients  “y.n.k  are  li sted  explicitly 
in  the  Appendix.  For  convenience  let  us  also  define 


a  *  —  0 

Y,n,k 


(5.14) 


for  k  <  0  or  k  >  2n. 


PROPOSITION'  (5.2': 


Let  0  <  n  <  N.  Then 


(5.17) 


t 


where  the 

(5.10). 


J 

+  y 

[2a.  fk<) 

+bi-l  (i-l)1  eY.n- 

Ai 

i 

i  > 

K 

+  l 

[a  ( 

i) 

+bi-l  (i-l)  +ci-2 

ik 

l  v 

upper 

summation 

limits  1,  J  and 

PROOF: 


up  to 


The  proof  is  identical  tc  that  of  Proposition  (5.1) 
the  point  where  in  place  of  (5.13)  we  would  obtain 


jJ=l  LkFy,m-k 


2m- 1 


6  H  „ 

p  y,m,n  y+m-2-n 


n 


for  m  >1  where  the  coefficients  6  are  given  by  (5.17). 

In  particular  for  m  =  1  Y,m,n 


eY,l,0  =  bo 


Because  of  the  different  boundary  and  initial  data  in  the  present 
case  however  we  now  use  Proposition  (4.1)  instead  of  Proposition 
(4.2)  and  obtain 


F  ( x » t )  =  -i— ^ 
Y.m  rr~  ' 


which  can  be  rewritten  as 


2m- 1 

—  Q 
m,n  aQ 

l 

n=0 

PY. 

2m-l 

y 

,m 

,  n+1 
1  X 

L 

n=0 

,n  ao  k=l 

a  as 

2 

2m 

X 

8  y  ,m,k' 

L 

k=l 

n+1  ,  n+2-k 


. (x,a  t) 


1  .  ,  (x,a  t) 

Y+m-k,k  o 


where  the  B  are  given  by  (5.16)  above  for  m  =  1,2,...,N. 

This  completes  tne  proof. 

For  n  =  0,1,2  the  coefficients  8  ^  are  given  explicitly 

in  the  Appendix.  Aswith  the  ay,n,k  use  define  for  conven- 

ience  B  .  =  0  (5.18) 

Y,n,k 
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for  k  <  0  and  k  >  2n. 

The  following  lemna  permits  us  to  characterize  the  functions 
directly  in  terms  of  the  functions  E  .  and  F  ,  . 


# 


E  n  ,  and  ,  directly 
y  ,  k  Y  *-k 


•»k 


.k 


LEMMA 


defined  by 


If 


{un:  n  =  0,1 . N}  E  0N(D,nN  )  and  uj  is 


u*(x,t)  =  (-1)°  u  *(x,t) 


then  {  u^:  n  =  0,1 . N  }  ef>N(  n  N )  i  here  u*(x,t)  =  un(-x,t) 


PROOF  : 


Let  us  transform  (5.2)  by  x  =  -y  and  =  -D  . 

Then  we  have 

lDraoDylun(-y*t)  =  Xlak  1Z£r"  Dy~bk-1  T^TTT  Dy 


k-2 


(  ) 

+ck-2  (k-2)!  1  un<‘>'>t 


“  .1.  *bk-l  !^td. 


k=l 


This  implies  that 


k-2 


+ck-2  '(£-’?)  !'',un(  y,t)‘ 


[D -L  ]u*  =  T  (-1)K  L.  u*  .. 

t  o  n  k  n-k 


Multiplying  across  by  (-1)  we  then  obtain 

n 


[D -L  ]u#  ~  . 

1  t  o J  n  k  n-k 


which  concludes  the  proof. 

Clearly  then  we  may  set 


and 


E"  =  (-1)  E*  „ 

Y,  n  Y,n 

F#  =  (-1  )nF*  n  =  0,1 . N 

Y  »n  Y,n 


(5.19) 

(5.20) 


since  these  functions  will  satisfy  the  correct  initial  and 
boundary  conditions. 
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The  coefficients  aYtn,k  and  BY,n,k  have  been 
listed  explicitly  in  Appendix  B  for  n  =  0,1, V.  ’because  they 
quickly  become  tedious  to  calculate  as  n  increases.  Perhaps 
symbol  manipulation  programs,  which  have  come  into  vogue, 
would  be  quite  useful  in  this  regard. 


R  h  M  A  R 1- 


Although  F  v.n  and_  F$n  were  originally  defined 
only  in  the  quarter  planes  Q  and  Q*  their  representations  (5.15) 
and  (5.20)  are  defined  throughout  fT.  This  permits  these  functions 
to  be  analytically  continued  (with  respect  to  x)  into  all  of 


H  in  such  a  manner  that  (5. 21,  remains  valid  in  the  extended 

domain  H.  Thus  {F  1  and  { F  *  }  are  of  class  PM(H,nM). 

lYtn  1  y*n  N  N 

By  setting  x  *  0  in  (5.7)  and  using  (3.16) 
and  (3*14)  we  have 


/rr  Y+n 

1  /aot 

E  (0,t)  = 

Y*n 

?  a  y  ,n,o  (  (y-hy)/2  )  ! 

(5.19) 

E^  (0,t)  = 
Y.n 

/ — -  Y+n 

,  ,,n  /at' 

(-1)  o 

i  a  y ,n,o  ( (y  +n)/2) ! 

0  in  (5.15) 

we  have 

1 

ro 

Y+n 

rY+niB  x  ■  i 

/a~  Y  k=0 
o 

Ik  J  B  y  .n,k  (Y+n) !  * 

(5.21) 


(5.21)* 


x  >  0  (5.22) 


„  „  i *n  2n  (-x)Y+n 

-  7=^  J0  (T)  rnsrr- 


x  <  0  (5.22)* 


R  !•!  M  ark 


For  any  a  >  0,  substituting  (3.22)  into  (5.7) 
and  (5.15)  easily  leads  to  the  following  identities 


EYtk(x,t)  =  /S  Y+k 


E  ,  (x/  /a7  t/a) 

T  t  K 


Fy  k(x,t)  =  /a~Y+k  F^  k  (x/ /a,  t/a) 


(5.23) 

(5.24) 
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VI.  FORMAL  EXPANSION  PROCEDURES— CONTINUED 

At  the  end  of  Section  11  we  obtained  a  system  of„  equations 
and  initial/boundary  data  for  the  individual  terms  in  the 

assumed  form  (2.9)  of  an  asymptotic  expansion  for  u(x,t). 
Comparing  equations  (5.2)  and  (5.4)  -  (5.6)  which  defined 
E^k  and  Fv  ^  in  Section  V  with  the  equations  (2.11)  and 
( 2  .*13 )  -  (2.15)’  formulated  for  we  immediately  obtain  the 

solutions 

Va.T)  =  EY,k{°*T) 

for  problem  ( I  VP )  ^  and 

uk(0. t)  =  Fy  k<o , t) 

for  problem  (BVP)Q  where  Ey  ^  and  Fy  ^  are  defined  by  (5.7) 
and  (5.15)  respectively.  'Substitution*  back  into  the  formal 
asymptotic  expansion  then  yields 

oo 

u(0,  t)  ~  eY  l  c 
k=0 

and 

00 

u(o  . T  )  ~  eY  J  e 

k=0 


By  re-expressing  these  expansions  in  terms  of  the  original 
variables  x  and  t,  related  to  a  and  T  by  (2.7) ,  and  truncating 
to  N  terms  we  obtain  the  following  functions 


VN  (x,t)  = 

N 

l 

eY+*<  E  .  ( (x-x  )/e  , 

t/  e  5  ) 

Y’xo 

k=0 

Y »  K  O 

WN  (x,t)  = 

N 

l 

F  .((x-x  )/e. 

t/e  2  ) 

Y’xo 

k=0 

y  *  k  0 

Recalling  identities 

(5.23),  and  (5.24) 

these  can 

be  written 

more  simply  as 

VN  (x,t)  = 

N 

l 

E  ,  (x-x  .t) 

(6.1) 

Y’Xo 

k=0 

Y,k  o 

E  .  (a  ,  t) 

Y»k 


Fr,k(0’T)* 
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(6.2) 


WN  'x.t) 
Y*x 


At 


-  y 


N 


o  k=0 

or,  in  the  more  expanded  form. 


l  FY,k(x_x«ft) 


M  N  2k 

Vy  „  »  l  l  a  k  i  H  +k  i  i  (x"xo’aot)  (6*3) 

Y.x0  k=0  j=0  Y»X,J  Y+x  J » J  0  ° 

M  _  N  2k 

WN  ( x , t )  =  2  /  a  At  Y  l  l  6  k  -H  (x-x  , a  t)  (6.4) 

Y»x0  0  k=0  j=0  Y»X»J  Y4k-J»l  o  o 

For  k  =  0,1,2  the  coefficients  k  i  anc*  6  y.k.i  are  g^ven 
explicitly  in  the  Appendix;  more  'generally  they’  can  be  deter- 


and 


Consideration  of 

the  converse  problems 

(IV?)  ff  and 

lead  to  the 

analogous 

5  expansions 

o 

V*N  (x,t) 

II 

.  (x-x  ,t) 

(6.5) 

Y*xo 

k=0 

Y,k  o 

W*N  (x.t) 

=  ATY 

I  v(x-x  ,t) 

(6.6) 

Y,xo 

k=0  Y,k 

which  can  also  be  written,  using  (5.19)  and  (5.20),  as 


#N  (x,t) 

N  L. 

=  l  (-l)k 

2k 

l  « 

Y*xo 

k=0 

j=0  Y’ 

^  (x.t) 

=  2/a  At  “ 

N 

Y  l  (-1) 

Y*xo 

o 

k=0 

(6.7) 


B  v.  ,H*  ,  .  (x-x  a  t)  (6.8) 

0  Y»k,J  Y+k-)*l  °*  ° 

The  expansions  V  f v  and  V^Xq  will  be  referred  to  as 
the  N-term  special  expansions  at  xr  with  respect  to  the  operator 
L  or,  more  briefly,  the  special  expansions.  The  expansions 

WN  and 

vxo  1-  •xo 

will  be  termed  the  N-term  right  and  left  hand  boundary  expan¬ 
sions  (respectively)  at  x  with  respect  to  the  operator  L7 

The  values  taken  on  by  the  special  expansions  at  x=Xo  are  of 
special  importance  and  can  be  obtained  by  substituting  (3.16) 
into  (6.3)  and  (6.7)  with  x  =  xQ.  We  have 
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and 


UN  ,  .v  1 

V  ( x,  t )  =  -K 

Y,xo  o 


v  (x  ,t)  =  l 
Y,xo  o 


N 

l  “y.k.o 
k=0 


N 

I 

k=0 


1  (-D  a 


/a”T  Y+k 
( ( Y+k  )'72")"T 


’’a  t  Y 
o 


+k 


Y,k,o  { (Y  +k)/5) ! 


(6.9) 


(6.9)# 


In  the  next  two  sections  we  shall  develop  uniformly  accurate 
approximations  for  the  solution  of  the  Cauchy  problem  and 
the  first  initial-boundary  value  problem  for  equation  (2.1) 
using  the  expansions 


Y.X, 


W 


N 


rm 


and  W 


Y  » Xr 


Note  that  these  functions  are  well  defined  whenever  the  Taylor 
coefficients 


ak 

k  =  0,1,.. 

.,N 

bk 

k  =  0,1,.. 

.  ,N-1 

ck 

k  =  0,1,.. 

. ,  N-2 

exist  at 

x  =  x0  .  In 

order  to  establish  rigorous  error  bounds 

for  the 

approximations 

of  the  next 

two  sections,  however,  it 

will  be 

necessary  to 

make  slightly 

stronger  assumptions  in 

the  form 

a(x) 

e  CN+a(I) 

(A) 

b(x) 

e  CN+1+°(1) 

(if  N  >  1) 

c(x) 

e  CN+24°(1) 

(if  N  >  2) 

where  N 

>  0  is  an  integer,  0  £  a  £ 

1  and  1  is  an  interval 

with  x0  e  1;  if  xo  corresponds  to  a  boundary  of  the  solution 
domain  for  (2.1)  then  x  Q  is  taken  as  an  end  point  of  1,  other¬ 
wise  x  0  must  be  an  interior  point  of  1. 


31 


VII.  THE  CAUCHY  PROBLEM 


The  Cauchy  problem  is  formulated  by  equations  (2.1)  and 

(2.2) 


Dt  u  =  L  u 


x  e  R,  t  >  0 


u(x,o)  =  f  (x) 


X  e  R 


If  a,  b  and  c  are  bounded  and  Holder  continuous  and  f  is 
continuous  then  this  problem  is  known3  to  have  a  unique  solution 
in  the  class  of  functions  satisfying  (2.6).  If,  in  addition, 
a,  b  and  c  are  sufficiently  differentiable  then  the  expansions 
developed  in  the  previous  sections  can  be  used  to  approximate 
u(x,t)  with  uniform  accuracy.  A  rigorous  analysis  of  this 
problem  was  undertaken  in  Reference  1  and  in  this  section  we 
wish  to  present  the  main  results  in  a  convenient,  abbreviated 
form  without  proofs. 


If  f ( x )  is  sufficiently  smooth  then  only  the  regular  expansion 
(2.3)  is  needed  to  approximate  u(x,t),  as  we  see  in  the  following 
two  theorems.  The  first  is  a  global  result  requiring  stronger 
hypotheses  than  the  second  which  is  a  local  result.  In  the 
following  we  shall  use  the  notation 


N  =  [N/2]  =  greatest  integer  <  N  /2 

E_c»- 

(N+l ) 


Let  a,  b,  c  e  CN+Q(R)  and  f  e  CN+2+“(R).  If 


L'"Ti'  f(x) 

is  bounded  for  all  x  e  R  then  there  exists  a  constant  K  >  0  and 
a  At  >  0  such  that  “ 


|  u ( x , t )  -  U^(x,t)  I  <  K  t^+1 


(7.1) 


uniformly  in  R  x  [o,At]. 


k. '  l  ■ 


See  Theorem  (2.2.1)  of  Reference  1. 


Let  l  C  R  be  an  open  interval,  let  (A  )  hold 
in  1  with  N  >  2,  a  >  0  and  let  f  e  C^+0I(I).  Then  there  exists 
a  At  >  0  anH  for  each  compact  set  SCI  there  exists  a  constant 
K  >  0  such  that 


|  u ( x , t )  -  U¥(x,t)  |  <  K 
uniformly  in  S  *  [o,  At]. 


(7.2) 


See  Theorem  (2.2.4)  and  Corollary  (2.3.1)  of  Reference 


In  many  practical  problems  the  initial  value  function  is 
not  smooth  but  can  be  accurately  represented  by  a  piecewise 
smooth  function.  Within  a  bounded  interval  1  this  can  then 
be  expressed  as  the  sum  of  a  smooth  function  and  (left  and 
right  hand)  jump  functions;  that  is 

K 

f(x)  =  f(x)  +  £  a ,  h  (x-x.)  (7.3) 

k=l  K  Yk  K 

K* 

+  I  h%  (x-xk> 

k=l  k  Yk  K 

where  7(x)  is  smooth  in  1,  x  ,  x*  el  and  a  ,  a*  t  0;  if 
no  jumps  occur  then  the  summation  terms  are  omitted. 

Since  our  problem  is  linear  and  since  the  two  previous 
Theorems  apply  to  f(x)  then  it  only  remains  to  treat  the  two 
"canonical"  cases 


f ( x )  =  h  (x-x  ) 

and  Y  ° 

f(x)  =  h*(x-x  ) 
Y  O 

We  shall  use  the  notation 


u 


and 


u 


*’xo  Y’xo 

respectively  to  denote  the  solutions  to  the  Cauchy  problem  with 
these  particular  choices  of  initial  data.  The  special  expansions 
(6.1)  and  (6.5)  were  constructed  in  the  last  four  sections  for 
the  express  purpose  of  approximating  these  functions;  the 
following  theorem  indicates  where  and  how  well  this  is  achieved. 


I  HI  >  >R| 


xQ  e 
exist 


Let  y  ^0,  let  1  C  R  be  an  open 
l  and  let  (A)  hold  in  l  with  N  >_  2,  a  >  0. 
a  constant  K  >  0  and  a  At  >  0  such  that 


interval 

Then 


with 

there 


u 


r.x, 


(x,t)  -  x  (x,t)  i  <  k  /rY+N+a 
>  Y’  o 


(7.4) 


uniformly  in  (-»,  xQ]  x  [o,  At]  and 

|u#  ( x « t )  —  V#Nv  (x,t)  |  <  K  /tY+N+“  (7.4)# 

*'xo  Y’x0 

uniformly  in  [xq,®]  *  [o,  At]. 


See  Theorem  (2.3.1)  of  Reference  1. 


This  result  is  not  altogether  satisfactory  since  it  still 
does  not  enable  us  to  write  an  approximation  for 


x  u  #  x 

y,xo  Y«Xo 

which  is  uniformly  valid  in  regions  such  as  S  x  [o,at]  where, 
S  is  a  compact  set  containing  x0  as  an  interior  point.  However, 
in  the  case  where  y  is  a  non-negative  integer,  say  y  =  n 
0,  we  can  circumvent  this  limitation  using  the  simple  identity 


(x-x  ) 
o 


h  (x-x  )  +  (-1)  h* ( x-x) 
no  no 


(7.5) 


This,  together  with  linearity  and  uniqueness  for  the  Cauchy 
problem,  implies  that 


+  (-l)n  u  f. 


(7.6) 


where  u  denotes  the  solution  having  the  polynomial  initial 

values  n,xo 


Un,xo(x*o)  =  U"Xo)  /n! 


Thus,  to  construct  an  expansion  for  x  valid  in  the  region 
x  >  xQ  we  need  only  re-arrange  (7.6)  in’tcP  the  form 

u  =  u  (~l)n 

n,xo  n’x0  n'xo 

The  first  term  on  the  right  can  be  approximated  using  its  N- 
term  regular  expansion 


l  L  x  [ (x-x  )n/n ! ]  tk/k! 
k=0  0 


(7.7) 


with  an  error  term  of  order  /tN+a  (Theorem  (7.2));  the  second 
term  can  be  approximated  for  x  xQ  using  the  N-n  term  special 

expansion  V  with  an  error  of  order  ^+a  (Theorem  (7.3)). 


We  are  thus  led  to  define  the  following  N-th  order  uniform 
expansion  for  u 


Similarly  we  can  define  the  N-th  order  uniform  expansion  for 

u  *  by  - - 

— n,  x  3 


y*N  (x.t) 
n-xo 


(-l)n  [\j*  (x.t)  -  V*)  "  (x.t)]  x<  X 


r-11  ln.: 

I  V*>.t 


n.x 


(7.8)# 


x  £  x 


We  can  now  state  the  following  result 

Let  (A)  hold  with  N  >  2  in  the  open  interval 
1  C  R  and  let  xq£  1.  Then  there  exists  a  At  >  0  and  for 
each  compact  subset  S  C  1  there  exists  a  constant  K  >  0  such 

that 


lu 


n.x. 


(x.t)  -  y  ^  (x.t)  |  <k  /t^ 


+  a 


n,  x 


and 


u*  (x.t)  -y*N  (x.t)  |  <  K  /IN+a 
n, x^  n ,  x  1  - 


(7.9) 


(7.9)# 


o  o 

uniformly  in  S  *  [o,  At]. 


Follows  directly  from  the  definition  of  the  uniform 
expansions  using  Theorem  (7.2)  and  (7.3).  This  theorem 

corresponds  to  Theorem  (2.4.1)  of  Reference  1. 


Unfortunately  we  have  not  been  able  to  construct  an  analogous 
uniformly  valid  expansions  for 

# 


u 


Y .  X 


and  u 


Y.x, 


O  '  o 

in  the  cases  where  y  is  not  an  integer. 


The  best  that  can 
dity  (x  £  xD  and 
x  i  xD)  for  the  special  expansions  can  be  enlarged  slightly 
to  include  the  parabolically  shaped  domain  (x-xD  )  2  <  4  Pt, 
0  £  t  <  i  t  where  P  >  0.  This  was  included  in  the  treatment 

of  the  Cauchy  problem  in  Reference  1. 


Although  the  restriction  to  integer  orders  does  reduce  the 
generality  of  our  theory  somewhat  it  is  still  sufficiently  flexible 
to  treat  most  physically  interesting  cases.  It  does  permit 
the  treatment  of  any  initial  value  function  expressed  as  a 
piecewise  polynomial,  for  instance.  Furthermore,  if  we  limit 
our  consideration  to  jumps  of  integer  order  then,  in  view  of 
identity  (7.5),  the  general  form  (7.3)  of  initial  value  functions 
can  be  rewritten  using  right  (or  left)  hand  jumps  only.  The 
following  theorem  thus  represents  our  most  general  result  for 
the  Cauchy  problem. 
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hold  in  1  with  N 
of  problem  (2.1), 
written  in  the  form 


Let  1  C  R  be  an  open  interval  and 
>^2,  a  >  0.  Let  u(x,t)  denote  the 
(2.2),  (2.6)  and  suppose  that  fix) 


let  (A) 
solution 
can  be 


K  N 

f(x)  =  Tix)  +  l  l  «kihi^x_xk) 
k=l  j=0  K]  ]  K 

where  f(x)  e  Cw  (1)  and  x,  e  1.  Then  there  exists  a  At  >  0 

and  for  each  compact  set  SCI  there  is  a  constant  K  0  such 
that 


mm  K  N 

|u(x,t)  -  T5N(x,t)  -  l  £  a  ( x , t )  |<  K  /tN+ “  (7.10) 

k=l  j=0  3,xk 

uniformly  in  S  *  [o,  At];  here  Y.  is  defined  by  (7.8)  and 
is  the  N-term  regular  expansion  for  initial  values  f(x): 

U^x.t)  =  I  L(k)  T(x)  tk/k! 
k=0 

with  N  =  [N/2] . 


VIII.  INITIAL-BOUNDARY  VALUE  PROBLEMS 

Consider  the  following  initial-boundary  value  problem  in 
the  domain  D  =  (x,,  x2)  x  (o.^t) 


Dtu  = 

Lu 

xo  <  X 

<  Xj,  O  <  t<4t 

u(x,o) 

=  f  ( x ) 

X„  <  X 

-  X1 

(BVP) 

o  — 

u(x,  ,t) 

=  g(t) 

o  <  t  < 

A  t 

u(x  2,t) 

=  h(t) 

o 

A 

r-f 

1  A 

A  t 

Using  linearity  this  can  be  broken  down  as  usual  into 
component  subproblems  by  treating  each  of  the  functions  fix), 
g(t)  and  h(t)  in  turn,  assuming  the  remaining  two  function 
to  vanish.  We  denote  the  problems  specified  in  this  manner 
by  (BVP)f  ,  (BVP)g  and  (BVP)^.  The  latter  two  problems 
are  nearly  identical,  of  course,  and  any  statement  concerning 
one  can  be  directly  transformed  into  an  analogous  statement 
concerning  the  other.  In  the  following  discussion  we  shall 
first  consider  problems  (BVP)  and  (BVP)h  and  then  turn  to 
problem  (BVP)f.  ®  ’’ 


$ 
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The  most  important  particular  forms  of  problems  (BVP) 
and  (BVP)  ^  are  those  with  the  canonical  boundary  data  ® 


g(t)  =  hy/2(t) 
for  problem  (BVP)^  and 
h(t)  =  hy/2 


o  <  t  <  A  t 


O  <  t  <  A  t 


(8.2) 


for  problem  (BVP)^.  The  boundary  expansions 


were  defined  by  (6.2)  and  (6.6)  with  precisely  these  problems  in 
mind  and  we  have  the  following  result. 

Let  y  ,  N,  a  >  0,  where  N  is  an  integer,  and 

let  (A)  hold  with  N  +  a  >  *0  in  an  interval  [xi  ,  Xj  +  d] 

for  some  d  >  0.  If  uq  denotes  the  solution  of  (BVP)  such 

that  f(x)  =  h(t)  =  0  and  (8.1)  holds  then  there  exist  a  At  >  0 

and  a  constant  K  >  0  such  that 


|  u  ( x ,  t )  -  W#N  (x,t)  I  <  K  /tY+N+“ 

8  t  iX, 


(8.3) 


uniformly  in  D.  Similarly  if  (A)  holds  with  N  +  a  >  0  in  [x2  -d, 
x2  ]  for  some  d  >  0  and  if  denotes  the  solution  of  (BVP) 

such  that  f ( x )  =  0,  g(t)  =  0  and  (8.2)  holds,  then 


I  u  (x,t)  -  w“  „  (x,t)  |  <  K  /?Y+N+a 


(8.4) 


uniformly  in  D. 


1BS2)  See  Theorem  (3.3.1)  of  Reference  1. 

The  intuitive  reasoning  behind  this  result  is,  first,  that 
by  their  very  construction  the  boundary  expansions  will 
accurately  approximate  Ug  and  near  the  respective  boundaries 
x  =  x  i  and  x  =  xa,  and  secondly,  that  these  expansions,  and 
these  derivatives,  decay  exponentially  away  from  the  first 
boundary  and  thereby  effectively  satisfy  the  homogeneous  equa¬ 
tion  (2.1)  and  vanish  at  the  opposite  boundary.  The  theorem 
is  established  by  incorporating  these  ideas  into  a  formal  appli¬ 
cation  of  the  maximum  principle  for  parabolic  differential  equa¬ 
tions. 

Theorem  (8.1)  can  clearly  be  generalized  to  include  functions 
g(t)  and  h(t)  which  are  finite  linear  combinations  of  functions 
h  Y /2  ( t ) .  Furthermore,  if  g(t)  (or  h ( t ) )  possesses  a  jump  dis¬ 
continuity  in  the  form 


V 


* 


0  0  <  t  <  t 

—  —  o 

/t=t^/(y/2)l  t  >  tQ 
where  0  <,  t  £  At  then  can  be  approximated  by  the  function 
j  0  0  <  t  <  tQ 

g  *  W^N  (x,t-t  )  t  <  t  <  At 

y ,  X  O  O  —  — 

Thus,  Theorem  (8.1)  can  be  directly  extended  to  include  any 

piecewise  continuous  boundary  values  of  the  form 

g(t)  K 

or  }  =  l  h  ,  (t-t.)  +  e(t)  (8.5) 

h(t)  k=l  K  Yk/2  k 


g(t)  = 


h  /0(t-t  ) 
y/2  o 


where  0  _<  y  <  N  and  e(t)  is  an  error  term  such  that 

|  e(t)  |  <  K  /tN+a 

for  some  constant  K  0.  In  this  more  general  case  however 
the  term  on  the  right  hand  sides  of  (8.3)  and  (8.4)  should 

be  replaced  by 

k  /rN  +  a 

since  the  lowest  possible  value  of  y  ,  y  -  0,  must  be  used. 
We  shall  not  formulate  an  additional  theorem  to  express  these 

results  since  they  are  straightforward  consequences  of  Theorem 
(8.1).  This  concludes  our  treatment  of  problems  (BVP)  and 

(BVP) ^  and  we  turn  now  to  problem  (BVP)f.  s 

The  regular  and  special  expansions  which  we  applied  to 
the  Cauchy  problem  in  Section  VII  are  still  meaningful  in  the 

context  of  problem  (BVP)  ,  provided  f(x)  and  the  coefficients 
a,  b  and  c  are  sufficiently  differentiable.  In  fact,  for  sufficient¬ 
ly  small  times,  these  expansions  will  accurately  approximate 
the  solution  of  (BVP)  except  in  regions  close  to  the  boundaries 
x  =  x  ,,  x2.  Theorems  (7.2)  and  (7.3)  of  the  last  section 
remain  valid  with  only  minor  modifications  and  may  be  stated 
as  follows: 


UBHilSHHiH  Let  I  C  (x,  ,  x  2  )  be  ari  open  interval,  let 
(A)  hold  in  I  with  N  +  a  >1  and  let  fe  C  ^(I).  If  Uf  denotes 
the  solution  of  problem  (BVP)  when  g(t)  =  h(t)  =  0  and  if 

TT 

lj  (x,t)  =  l  L(k)  f ( x )  tk/k!  (8.6) 

k=0 
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where  fT  =  [N/2]  then  there  exists  a  At  >  0  and  for  each  compact 
subset  S  C  I  there  exists  a  constant  K  such  that 

|  uf(x,t)  -  UN(x,t)  |  <  K  (8.7) 

uniformly  in  S  x  [o,At]. 
mm  See  Theorem  (3.2.4)  of  Reference  1. 


fill  -.1 K  t  \! 


Let  1  C  (x  i  ,  x  2  ) 


(A)  hold  in  1  with  N  +  a  >  0  and  let 
the  solution  of  (BVP)  such  that  g(t) 
that 


be  an  open  interval,  let 
x0  e  I.  Let  uf  denote 
=  h(t)  =  0  and  suppose 


f ( x )  =  h^(x  -  xQ) 

for  some  y  >_  0  and  xi  <  x  <  Xj  .  Ify  +  N+  o>l  then  there 
exist  a  constant  K  >  0  and  a  At  >  0  such  that 

| Uc( x , t )  —  VN  (x.t)  |  <  k/Fy+N+“  (8.8) 

1  Y  ,xo 

uniformly  in  [x,,  x0]  *  [o,At].  Similarly  if 


then 


uf(x,t)  - 


f  ( x ) 
,#N 


=  h*  (x-x  ) 
Y  o 


V""  ( x , t )  |  <  K  /t 

Y  .X 


Y  +n+  a 


uniformly  in  [x0,  x2l  x  [o,At].  Here  the  functions  VN 
vffN  are  defined  by  (6.1)  and  (6.5)  Y" 


y.x. 


See  Theorem  (3.2.5)  and  (3.3.3)  of  Reference  1. 


(8.9) 

and 


By  combining  Theorems  (8.2)  and  (8.3)  we  can  obtain  a 
result,  for  problem  (BVP)  f,  almost  identical  to  Theorem  (7.5) 
in  which  initial  values  of  the  form 

K  N 


f(x)  =  F(x)  +  l  l 

j«0 


akjhj(x-xk) 


Nfaksl 

are  treated,  where  f  e  C  (1)  with  1  C  (x,  ,  x2  ).  The  only 
modification  would  be  that  the  restrictions  on  N,  a  and  y  would 
now  be  N,  a,  y  >_  0  with  N  +  a  >  1.  We  omit  a  formal  statement 
of  this  Theorem  . 


Our  analysis  is  still  not  complete  however  since  we  have 
little  information  about  the  behavior  of  U  near  the  boundaries 
x  =  x  i  ,  x  2  ,  even  when  the  open  interval  in  Theorem  (8.2) 
and  (8.3)  can  be  taken  as  l  =  (x,  ,  x2  ).  The  dependence 
of  the  constant  K  on  the  compact  subset  S  C  l  is  a  critical 
weakness  since  it  means  that  no  single  constant  K  exists  such 
that  inequalities  (8.7)  and  (8.8)  hold  uniformly  in  all  of  D. 
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Let  us  consider,  therefore,  how  to  improve  our  approximations 
in  the  vicinity  of  x  =  x2. 

Since  the  closest  discontinuity  in  the  initial  values  f(x) 
(if  any  exists)  occurs  at  some  finite  distance  away  from  x 
=  x  2  then  we  may  suppose  that  f(x)  is  smooth,  say  of  class 
C  N+ a  some  neighborhood  of  x  =  x2  .  The  regular  expansion 
(8.6)  is  thus  well  defined  for  x  in  this  region  but  cannot 
approximate  Uf  accurately  since  along  x  =  x2  it  assumes  the 
(in  general)  non-vanishing  values 

U^(x2,t)  =  V  L(k)  f(x2)  tk/k! 
k=0 

In  effect  is  approximating  the  solution  of  the  wrong  problem, 
namely  one  for  which 

h(t)  =  l  L(k)  f(x2)  tk/k !  (8.10) 

k=o 


along  x  =  x2  instead  of  h(t)  =  0.  To  counteract  this  it  would 
seem  natural  to  adjust  U  "  by  subtracting  away  the  solution 
of  problem  (BVP)^  for  which  h(t)  is  given  by  (8.10).  But 
in  this  form  h(t)  is  simply  a  finite  power  series  in  t  so  the 
methods  already  developed  for  problem  (BVPk  apply  and  lead 
us  to  define  the  following  N-th  order  right  hand  boundary 
correction  for  problem  (BVP)^. 

wN(x,t)  =  T  l(k)  f(x2)  WN“2k  (x,t)  (8.11) 

k=0.  <Jk,x2 


Similarly,  we  define  the  N-th  order  left  hand  boundary  correc¬ 
tion  for  (BVP)f  by 


w*N(x,t) 


TT 

Jo 


(k) 


f(xt  )  W 


#N-2k 

2k,Xj 


(x,t) 


(8.11)# 


We  can  now  state  the  following  theorem: 


Let  N,  a  >0  with  N  + 


c  be  of  class  C  N+a 

N+  a 


in  (xi,  x2)  and  let  a,  b, 

(Xj  -  d,  x2]  for  some  d  >  0.  If  f  e  C 
if  Uj  denotes  the  solution  of  (BVP)  when  g(t) 
there  exist  a  At  >0  and  a  constant  K  >  0  such  that 


1,  let  (A)  hold 
[ x i ,  x,  +  d)and 
( [x,  ,  x2  ] )  and 
=  h(t)  =  0  then 


l(uf  -  + 


N 


w 


+  w 


#N 


)  ( x ,  t )  I  1  K  /t^ 
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(8.12) 


uniformly  in  D  =  [x^.x^]  *  [0,At]. 


See  Theorem  (3.3.6)  of  Reference  1. 


This  essentially  completes  our  discussion.  Uniformly  valid 
approximations  for  when 

f  ( x )  = 

or 

f  ( x )  = 

in  the  cases  where  y  =  n,  a  non-negative  integer,  can  be  obtain¬ 
ed  by  repeating  the  device  of  Section  VII  involving  identity 
(7.5)  and  applying  Theorems  (8.3)  and  (8.4). 


h  (x  -  x  ) 
y  o 

h*  (x  -  x  ) 

y  O 


Thus,  to  summarize  this  section,  we  can  construct  approxi¬ 
mate  solutions  to  problem  (BVP)  which,  for  sufficiently  small 
At  >  0,  are  accurate  to  order  At  +a  uniformly  in  D  whenever 
g(t)  and  h(t)  have  the  general  form  (8.5)  and 


f(x) 


f(x)  + 


N 


a  .  .h . (x 
kj  ] 


-  x, 


) 


where  f  e  C  +  a  ([x,,  x2]).  Each  of  the  terms  in  these  forms 
can  be  treated  separately  because  of  the  linearity  of  problem 
(BVP)  and  the  general  result  obtained  by  superposition. 
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APPENDIX  -  THE  COEFFICIENTS  «  „  .  AND  B  „ 

Y  § n  f k  y  a 

The  coefficients  aY,n,k  and  appearing  in  (5.7) 

and  (5.15)  respectively  can  be  obtained  from  formulae  (5.8) 
and  (5.16)  respectively.  Since  these  become  somewhat  laborious 
to  compute  explicitly  as  n  increases  we  include  here  the  results 
for  cases  n  =  0,1,2.  For  larger  n  it  would  seem  advisable 
to  utilize  a  symbol  manipulation  program  to  double  check  the 
results. 

In  the  following  formulae  let 

=  ai/ao  i  =  0,1,2 

Bi  =  bi/aQ  i  =  °,1 

*i  =  ci/ao  1  =  0 

n  =  0 


aY,2,0  y(  Y+2)  (2Y~1  )S>2'*' "J  Y  (y+2)Bj  +  -^y(y+2)C0 

+  Y2  (y2-4)S  J  +  Jy*Y-1)(y+2)s^0  +  i  y(y+2)\2 

\,2.l  -  -  bVz  V  W  si  <2^ira?  - 


_ a... 


■  i 


av92,2  =  zVl?  1  "  T5  (2y2+1)si2  -  ZY*lbo  +  Z  ^0 


Y» 
aY,2,3 


3-2 

a  Y.2,4  =  Z  ®L 

6y,2,0  =  0 


1  K2  +  I  Slbo  +  I  312 


„  1*  1  r  1-  3  ,2  1-r  It  2 

6  y,2, 1  =  ‘  S  S2  +  Z  bl  "  7  Co  +  3T  ai  ‘  Z  Slbo  +  -Bbo 


1 


6  y,2,2  =  i  S2  ’  \  *1  “  IB  al  +  Z  5lbo  +  Z  bo 


3  —2  .  1  —  t  ,lr  2 


Y. 


1  -  .3-2.3,* 


2,3  =  ”  I  a2  +  B  *1  +  Z 


1  o 


o  _  3-2 

By,2,4  "  Z  al 


For  completeness  we 

a  _  1,  and  3  ., 

T  ,ntK 

and  (5.17)  respectively: 


also  list  the  intermediate  coefficients 
n  =  1,2,  appearing  in  equations  (5.9) 


a  y  »1*0  bo 
“y.1.1  =  al 
3y,1,0  =  bo 
*  Y.1.1  =  ai 

rl  2  _  -r  l  t-2 

°y,2,0  =  ao[Z  Y  albo  +  +  CoJ 

\,2,1  •  aolZ  (^-1)Sl2  -  \  bo  +  Z  <2^3)alVbl> 

*y,2,2  =  V  _  \  312  “  \  albo  +  *2] 

3  y.2,3  =  aol  -  \  V1 
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